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Abstract 

We present theoretical predictions for a few phenomenologically interesting 
distributions in the semileptonic b — > u decays which are affected by Fermi 
motion. The perturbative effects are incorporated at the onedoop level and 
appear to be very moderate. Our treatment of Fermi motion is based directly 
on QCD, being encoded in the universal distribution function F(x). The 
decay distributions in the charged lepton energy, invariant mass of hadrons, 
hadron energy, and q 2 are given. We note that typically about 90% of all 
decay events are expected to have Mx < Md; this feature can be exploited 
in determination of IV^&I. 
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1 Introduction 



The measurement of heavy quark decay distributions is important in extracting key 
parameters of the standard model. Towards this end, theorists are faced with the 
tough task of separating strong binding effects, arising from strong interactions at 
large distances, from the relatively simple quark-lepton Lagrangian, known at short 
distances. Both perturbative corrections and little understood nonperturbative ef- 
fects are very moderate in beauty decays if one is interested in sufficiently inclusive, 
integrated characteristics [I]> 0- In experimental studies of b — ► utv decays, attempts 
to separate b — > c transitions without explicitly identifying charmed particles, often 
suggests using a narrow slice of the overall decay kinematics reflecting the low in- 
variant mass of final-state hadronic system accessible only in b — > u decays. The 
strong interaction effects are magnified in this case. This fact was realized already 
in early papers where the QCD-based treatment of inclusive decay rates was elabo- 
rated. In particular, it was pointed out in |2| that the QCD-based OPE itself leads 
to the emergence of the phenomenon of "Fermi motion" of the heavy quark inside 
a hadron, an effect introduced ad hoc much earlier || for the description of heavy 
flavor decays in phenomenological models. For example, the impact of Fermi motion 
in simple quark models on the invariant mass of hadrons, Mx, was also addressed 
in ||. The concrete manifestation of Fermi motion in QCD has some peculiarities 
P, H making it somewhat different from the simple-minded Fermi motion of quark 
models. The QCD treatment was later described in detail in a number of publi- 
cations |5], |6], |7|, |8|, ||. Conceptually, Fermi motion is similar to the leading- twist 
nonperturbative effects in DIS. Further details and more extensive references can be 
found in ||. 

In the present work we apply the methods of to a few experimentally impor- 
tant distributions in the b — > u semileptonic decays, namely, the distribution of the 
electron energy dT/dEi, the invariant hadron mass dT/dMx and the hadron energy 
dT/dEh] we also give a plot of dT/dq 2 . The first three distributions are directly af- 
fected by the Fermi motion. Our treatment is in certain aspects simplified compared 
to H; we will mention these elements and justification later. 

2 General Strategy — Including Fermi Motion in 
Heavy Quark Decay Distributions 

Strong interactions have two faces in 6-decays: one, associated with the hard modes 
of the gluon fields, is to produce perturbative corrections, the other, associated 
with the soft modes is to produce nonperturbative corrections, in particular Fermi 
motion. The way to separate the hard and soft modes is dictated by Wilson's OPE: 
the introduction of the hard separation scale \i. The 'identity' of the two regimes 
has no absolute sense and depends on the choice of \x. 

The hadronic part of the semileptonic decays is described by five structure func- 
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tions u>i(qo, q 2 ), only three of which are relevant for I = e or /i (for definitions, see 
flOf ). The effect of Fermi motion is encoded in the heavy quark distribution function 
F(x); note that the primordial momentum is normalized to A = Mb — so that 
we deal with the dimensionless parameter x. Then the expression for the structure 
functions is written in the following way: 

W (q , q 2 ) 1Caam ^- / [ q - ^qf^f, q 2 ) F(x) ( 1 - ^^==j dx ' 



(1) 

Here w pert is a parton structure function dressed with short-distance (perturbative) 
corrections. This expression coincides with those of Ref. || through leading-twist 
terms, summation of which yields the effect of the primordial momentum distribution 
of the heavy quark. The form of Eq. ([!]) is convenient since it closely follows the 
simple picture of the decay of the heavy quark boosted to the velocity xA/M B 
along the direction of q, with F(x) determining the probability of the corresponding 
initial-state configuration. 

The moments aj of F(x) are given by the expectation values of local heavy quark 
operators. In the adopted normalization 

a = 1 , ai = , a 2 = -=2 , a 3 = — =3 (2) 

3A 3A 

with 

A = ^W(^) 26 l 5 > » Pd = -^(B\bD,G^b\B) . (3) 

One can chose a particular functional form for F(x) and adjust a few parameters to 
fit the phenomenologically deduced moments. This procedure was undertaken in 
where the following ansatz was suggested: 

F{x) = 0(1 -x) e cx {l -x) a [a + b{l - x) k ] . (4) 

To construct the distributions we are interested in, we thus use the perturbatively 
corrected parton distributions w pert (go, q 2 ), add the effect of the nonperturbative 
distribution via the convolution of Eq. (0), and calculate the decay distributions 
with the full w(qo, q 2 )- This procedure, therefore, amounts to averaging the pertur- 
bative decay distributions over nonrelativistic primordial motion of the heavy quark 
governed by the distribution function F. 

This approach is simplified in one important aspect. In QCD, contrary to non- 
relativistic models, F(x) depends essentially on q 2 and on the final state quark 
mass even when perturbative effects are switched off. In particular, F(x) com- 
pletely changes when m& — \fq 1 ^ M x where M x is the invariant mass of the final 
hadronic system. For example, relations (0) are modified. In the presence of gluon 
bremsstrahlung Mx can become large, 3> Aqcd, even for b — ► u transitions. Certain 
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relations between the moments of F at different q 2 still hold,[] and we expect that 
the effects of changing F(x) are insignificant in what concerns single-differential 
distributions. 

Another complication of the QCD description is that F(x) together with A and 
local operators defining moments a.; are normalization-point dependent when per- 
turbative corrections to parton distributions appear. The perturbative structure 
functions are /i-dependent as well; only the convolution (|l|) is ^-independent. The 
most noticeable effect that arises is due to the running mass m&(//), and related 
variation of A(/x) @. 



3 Fermi motion in dT/dEi, dY/dMx, and dY/dE^ 

With the structure functions (P we arrive at the following lepton spectrum in b 
u£u decays: 



dT 
~dE p 



— - — E t - —^x Fix) 1 - — — dx 
dE I M B v ; I M B 



(5) 



For the two other distributions, dT/dMx and dT/dEh, we need to recall the kine- 
matic definitions: 



M\ = Ml + q 2 -2M B q = \m\ + q 2 - 2m b q ] + 2(m b - q )A + A (6) 



and 

From this we get 
dT 



E h = M B -q = m b -q + A 



(7) 



dM 2 x 2Mb 
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Aq 2 M 2 B 



-1/2N 
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<r 2 P"- rt I ' M 2 B + q 2 -M 2 x xA 
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and 



dT 
dE~h 



J dq 2 J dx F(x) ( 1 - ^ (l - q 2 /(M B - E, 



.fr 1/2 i x 



dq dq 2 



M, 



B 



(9) 



1 The concrete form of Eq. (|l|) is chosen to preserve the lowest moments. 
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4 Perturbative corrections 



Perturbative corrections are parametrically enhanced in the kinematics close to the 
free-quark decay. In particular, in b — > u (b — > s) transitions double log effects ap- 
pear. These enhanced corrections must be summed up in the decay b — > s + 7 where 
q 2 = 0, which was done in |J, including the effect of running a s . In semileptonic 
transitions the bremsstrahlung effects are much softer since the typical configura- 
tion corresponds to a significant invariant mass of the lepton pair. Therefore, we use 



instead the exact one-loop radiative corrections calculated in Ref. ||11|| ; the running 
of a s is thus discarded as well. It seems justified since even in the case of b — > s + 7 
these effects - although very important in purely perturbative calculations including 
integration over low-momentum gluons, - are marginally seen when one proceeds 
with the Wilson's OPE where the evolution of the coefficient functions and effective 
low-energy parameters is done on the same footing. We will demonstrate also that 
the results are fairly insensitive to variation of a s , which is another justification. 

On the other hand, limiting ourselves to the one-loop fixed-a s calculations, and 
keeping in mind the extremely weak numerical //-dependence of the obtained physical 



spectrum^] we formally put \i — which allows using the expressions of Ref. (TT| 
literally. This, however, necessitates using the perturbative one-loop value of the 
6-quark pole mass 

rh b ~ m b (fi) + — — \i , (10) 
y 7r 

likewise a similar subtracted value of fi 2 

frl - t^M ~ i—V 2 , (11) 

O 7T 



etc. (for more details, see review [0). We use a s /n ~ 0.1 in the analysis; at 
li w 1 GeV it corresponds to rhb ~ 4.82 GeV, and jl 2 , being about 0.1 GeV 2 smaller 
than fi 2 normalized at the scale 0.5-1 GeV. 

Although the perturbative distributions are known to one loop only, one can 
easily write them in exponentiated form (see, e.g. 0). This is technically convenient 
since the elastic peak affected by virtual corrections disappears and one can then 
deal with smooth perturbative distributions. One relatively simple form of the 
exponentiated corrections for the double distribution was suggested in Eq. (17) of 
Ref. |13| . To simplify the calculations we used this form of radiative corrections 
for dT/dMx and dT/dE h where double-differential distributions must be used. The 
numerical difference of exponentiation is negligible after incorporating the effect of 
the primordial Fermi motion. 



2 The residual //-dependence always remains in theoretical calculations due to their approximate 
nature. 
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5 Results and Discussion 



The results of calculation of the distributions are shown in Figs. 1-3 (dT/dEi), Fig. 4 
(dT/dMx) and Fig. 5 (dT/dE h ). Our basic set of parameters is m b = 4.82 GeV, 
p,^ = 0.2, 0.4 and 0.6 GeV 2 which correspond to 

[a, a, b, c] = [1.5, 3, 16, -3.3] , [0.5, 2, 0.76, -1.75] and [0.1, 1.2, 0, -1.1] , (12) 

respectively, in the ansatz for F(x) (we use k = 1); the value of a s is set to 0.3. 

To illustrate dependence on m b we show also the plot for dT/dE^ for rh b = 
4.73 GeV (Fig. 2). Figure 3 shows the dependence of the lepton spectrum on a s ; it 
is clearly rather weak. 

As in the case of b — > s + 7, the effect of the Fermi motion - where present - 
is more pronounced than perturbative modifications. As expected, it shows up in a 
much softer way in dT/dE^ compared to dT/dE^ in b — > s + 7. The situation with 
b — > s + 7 is in closer analogy with dT/dM x , but even here the effect is suppressed 
due to significant average invariant mass of the lepton pair, i.e. effectively smaller 
energy release and recoil momentum. 

Apart from smearing by Fermi motion of non-smooth parts of the partonic dis- 
tributions, the most significant nonperturbative effect in dT/dMx and dT/dEh is 
associated with the difference A = M# — m b in the value of Eh = Mb — qo- This 



nonperturbative effect was first pointed out and analyzed in [H . 

The predicted shape of the distributions depends to some extent on the adopted 
form of the distribution function F(x). We believe, nonetheless, that possible vari- 
ations are not significant as long as the same underlying parameters, m fe (/i) and 
/i 2 (/i) are employed. 

A more significant effect is possible from the next-to-leading twist operators, 
since the effective energy release is not very large. An example of such effects is the 
chromomagnetic interaction whose effect can be estimated pi, |9|] , for instance, as due 
to the final-state mass difference between 7r and p, or due to the initial-state mass 
difference for weak decays of B and B* - these effects are neglected in the analysis. 
They can cause certain shifts of the distributions which are predicted if smearing 
over the corresponding interval is done. For example, in the distribution over M x 
the necessary interval of smearing over M x constitutes, probably, about ±0.2 GeV 2 . 
In the regular parts of the distributions such a smearing is superfluous, of course. 

For the same reason the very beginning of the distribution over M x and Eh 
cannot be taken literally; in particular, one cannot deduce from it the exclusive 
decay rate into 7r or p, or any particular resonance. In the limit m b — > 00 the 
distance between the successive resonances in M x would disappear in the scale of 
(Mx); for the actual case of B it still produces a relatively coarse 'grid'. 

A related peculiarity of the presented plots is that they show a nonvanishing rate 
for M x or Eh below pion mass, which, from the OPE viewpoint, is completely an 
artefact of neglecting higher-twist effects. To get rid of this unphysical feature one 
can, for example, consider ad hoc the given distribution as the one over M\ — M 2 
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rather than over M\ - distributions modified in such a way are formally equivalent 
to the order in l/ra& expansion we work. 

The decay distribution dT/dq 2 is not essentially affected by strong interactions. 
At maximal q 2 particularly close to ml, nevertheless, these effects blow up. It is 
important to realize that it is not literally the effect of Fermi motion that shapes 
dT/dq 2 . For example, the integral over the large-g 2 domain is governed by the flavor- 
dependent expectation value of the six-quark operator and thus can be completely 
different in decays of and B°. The effects originating at q 2 —>■ ml were dis- 



cussed in detail in ||14|| . Here we remind that studying the difference of the decay 
distributions for and B° near maximal q 2 , at small Mx or Eh - or merely in 
the end-point electron spectrum - is a way to directly measure the four-fermion 
expectation values which are important for B physics. 

The decay distribution dT/dq 2 is shown in Fig. 6; except near the kinematic 
boundary, the deviation from the tree-level parton shape is very small. We ob- 
tained this distribution merely adding known nonperturbative 1/m 2 effects to 



the parton expressions calculated through order a s fllTl . 

An interesting observation we infer from the plots is that a significant fraction, 
~ 90%, of the decay events is expected to have M x < 1.87 GeV, i.e. lie below the 
charmed states. It is in contrast with the case of semileptonic spectrum where only a 
small fraction of the decays proceed to the domain above the kinematical bound for 
b — > c transitions, and even small measurement errors lead to sizeable bias. For the 
distribution in invariant mass, with a lower cutoff on Mx between 1.5 and 1.6 GeV 
(to allow for a possible leaking of higher-mass states due to experimental uncertain- 
ties), the majority of decays appears in the low-Mx region (80% for M x < 1.6 GeV 
and about 75% for M x < 1.5 GeV), and can be reliably calculated theoretically. On 
the other hand, this can possibly be determined in experiment. This would suggest 
a way for a trustworthy determination of \V U b\ with a relatively good accuracy. 

If a measurement of the distribution over M\ is possible, it yields the possibility 
to determine independently the hadronic parameters. In analogy to the b — > s + 7 
decays, the center of gravity of the distribution is sensitive to the b quark mass and 
its width determines fi 2 . Figure 4 shows a reasonable sensitivity of the distribution 
over M\ to nl. 



6 On a Possible Improvement of Theoretical Pre- 
dictions 

The presented analysis of the decay distributions is in many aspects simplified. 
While the b quark mass is large enough for fully integrated characteristics, it is not 
the case when the detailed differential predictions are attempted. It is not therefore 
clear that further refinements based on more advanced treatment of perturbative 
corrections and nonperturbative effects can yield an essential and trustworthy im- 
provement. The first candidate for the refinement is, obviously, inclusion of known 
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1/m 2 corrections which are not incorporated in the leading-twist effects summed up 
by the Fermi motion. The anticipated scale of these effects was mentioned above. 

Another direction for refinements is improving the perturbative description. The 
impact of the perturbative corrections on the distributions even in B decays appears 
to be smaller than those of nonperturbative effects. In particular, we expect a small 
effect from inclusion of higher-order a s -corrections. This applies, however, only 
when a proper treatment of the infrared domain is done - otherwise the running 
of a s , typically the dominant effect among the second-order corrections, apparently 
generates effects which blow up. This does not happen when one literally follows 
Wilson's procedure of constructing the OPE. Experience shows that this procedure 
is necessary already when the second-order corrections in b decays are addressed - the 
procedure enforces the safeguard from inconvenience of the miracleous compensation 



of significant effects coming from different sources [15], [L6[ . 

While introduction of the IR cutoff in calculating the perturbative coefficient 
functions is more or less straightforward in usual perturbative calculations, this is not 
so simple when exponentiation of soft/collinear effects is employed. A method ap- 
plicable to these problems was described in ; it combines a few essential elements: 
reproducing exact one-loop and (all-order) BLM-improved answer and, simultane- 
ously, proper exponentiation of singular double-log corrections. This feasible way is 
expected to yield more than enough accuracy in evaluating perturbative corrections. 

While such a treatment is necessary when one intends to determine underly- 
ing parameters of the heavy quark expansion like rrib and /i^ with the accuracy 
when their scale-dependence is essential, we anticipate a small overall impact of the 
higher-order perturbative corrections on various hadronic characteristics involved, 
for example, in measuring \ V u b\, well below the effect of the leading nonperturbative 
phenomena. 



7 Outlook 

One of the most important practical applications of the analysis of the b — > ulv 
distributions is the extraction of \V u b\. As was realized long ago jflTfl ? the theoretically 



cleanest way to determine it is from the total b — > utv width. In particular |15| 



, , /BR (B -> XJu )\ 2 /1.55ps\ 2 . , 

\V ub \ ~ 0.00415 i ^ 1 —\ 13 

1 1 V 0.0016 J V r B J K ' 

where theoretical uncertainties lie well below the level which is experimentally rele- 
vant now and in the near future. However, a completely model-independent exper- 
imental measurement of such a width is embryonic at present. 

Another extreme way, suggested long ago, was to consider decay events with 
Ei > {M B — M^)/2Mb = 2.31 GeV, where charm decays cannot contribute. The 
decay rate in this too narrow slice of kinematics, on the other hand, depends on 
poorly known details of strong dynamics. Within our calculation this fraction of all 
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decay events does not vary too significantly; nevertheless, more work is needed to 
put this observation on a more firmly grounded quantitative basis. 

There are many intermediate options whose advantages mainly depend on exper- 
iment. The larger the kinematic domain, the more it encompasses the parton- level 
kinematics, the better, in general, is the theoretical description of such an inclusive 
width. We found that the distribution over Mx studied in experiment is rather 
promising in this respect. A similar suggestion was made to look at the distri- 
bution over E h [|13|]. It is clear that an accurate determination of \V u b\ requires a 



detailed analysis of the kinematics of the decay events, which would ultimately lead 
to determination of just the double-differential distributions over q Q and q 2 . Which 
particular integral appears to be most suitable for extracting \V u b\ will eventually 
depend on the available experimental technique. 
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Figure 1: The electron energy distribution dT s \(b — > u)/dEi, arbitrary units. Long- 
dashed, solid, and short-dashed lines correspond to = 0.2, 0.4 and 0.6 GeV 2 . The 
6-quark mass m b = 4.82 GeV, a s = 0.3. 
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Figure 3: Effect of radiative corrections on dT/dEf. solid line shows a s = 0.3 and 
dashed line is a s = 0. Here fn b = 4.82 GeV and fil = 0.4 GeV 2 . 
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Figure 4: The invariant mass distribution dT/dM x ; long-dashed, solid, and short- 
dashed lines correspond to fx^. = 0.2, 0.4 and 0.6 GeV 2 . The 6-quark mass fhb = 
4.82 GeV, a s = 0.3. All distributions are normalized to the same total width. 
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Figure 5: The hadronic energy distribution dT/dE h in the same setting as in Figs. 1 
and 4. 
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